Low-temperature and high-temperature approximations for penetrable-sphere fluids. 
Comparison with Monte Carlo simulations and integral equation theories 
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The two-body interaction in dilute solutions of polymer chains in good solvents can be modeled 
by means of effective bounded potentials, the simplest of which being that of penetrable spheres 
(PSs). In this paper we construct two simple analytical theories for the structural properties of 
PS fluids: a low-temperature (LT) approximation, that can be seen as an extension to PSs of the 
well-known solution of the Percus-Yevick (PY) equation for hard spheres, and a high-temperature 
(HT) approximation based on the exact asymptotic behavior in the limit of infinite temperature. 
Monte Carlo simulations for a wide range of temperatures and densities are performed to assess the 
validity of both theories. It is found that, despite their simplicity, the HT and LT approximations 
exhibit a fair agreement with the simulation data within their respective domains of applicability, 
so that they complement each other. A comparison with numerical solutions of the PY and the 
hypernetted-chain approximations is also carried out, the latter showing a very good performance, 
except inside the core at low temperatures. 



PACS numbers: 61.20.Gy, 61.20.Ne, 61.20.Ja, 05.20.Jj 



I. INTRODUCTION 

Since Wertheim and Thiele's analytical solution 
of the Percus-Yevick (PY) integral equation for hard 
spheres (HSs), a great deal of effort has been devoted 
to (approximately) describe the equilibrium behavior of 
simple liquids. Potentials capturing short-range harsh re- 
pulsion (with or without an attractive contribution) have 
deserved particular interest. Recently, the mainstream 
of attention has been transferred to more complex liq- 
uids, such as colloidal dispersions. Interestingly enough, 
coarse-graining techniques reveal that apparently toy in- 
teraction models, some of them not diverging at zero sep- 
aration, can serve as excellent effective potentials in soft 
condensed matter. 

While the effective interaction between two sterically 
stabilized colloidal particles can be accurately modeled 
by the HS potential [2], the effective two-body interac- 
tion in other colloidal systems can be much softer. For in- 
stance, the interaction potential for star polymers in good 
solvents diverges only logarithmically for short distances 
[3, Q . In the case of dilute solutions of polymer chains in 
good solvents, the distance between the centers of mass 
of two chains can be smaller than the sum of their respec- 
tive radii of gyration 0], so that the effective two-body 
potential is bounded. This implies the possibility of par- 



ticles to form clusters, which makes the problem greatly 
non-trivial. A simple class of bounded potentials is de- 
scribed by the purely repulsive generalized exponential 
model [1, S @ parameterized by an index n. In the par- 
ticular case n = 2 one recovers the Gaussian core model, 
which has been widely studied 0, i, i, [loLEl, El, [H, El ■ 
This model exhibits re-entrant melting [llj and the same 
is expected if n < 2 [1, Q. On the other hand, a clus- 
tering transition, which has been analyzed in detail for 
n = 4 0, i], is predicted ff n > 2 0, i, EH- In the hmit 
ri — i- oo, the generalized exponential model reduces to 
the penetrable-sphere (PS) model 



r < a, 
r > a. 



(1.1) 



This model has been extensive! 
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25|, |2a 123, 123, |29|, |30| . Density-functional theory 
predicts a freezing transition to fee solid phases 
with multiply occupied lattice sites. The existence of 
clusters of overlapped particles (or "clumps") in the PS 
crystal and glass was already pointed out by Klein et al. 
jl§ |. who also performed Monte Carlo (MC) simulations 
on the system. In the fluid phase, the standard inte- 
gral equation theories are not reliable in describing the 
structure of the PS fluid inside the core at low temper- 
atures [l^ [2lj : the number of overlapped pairs is over- 
estimated by the hypernetted-chain (HNC) theory, while 
it is strongly underestimated by the PY theory. Other 
more sophisticated closures l2ll.l25ll , as well as Rosenfeld's 
fundamental-measure theory [22| , are able to predict the 
correlations functions with a higher precision. 

In any case, the classical theories (PY and HNC), as 
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well as the alternative ones (2l|, [2^, l25| , require numerical 
work to get the correlation functions, such as the radial 
distribution function g{r) or, equivalently, the cavity (or 
background) function y{r) = e'^^'^^^''^'^ g{r), where fcs is 
the Boltzniann constant and T is the temperature. 

In general, the functions g{r) and y{r) convey the same 
physical information. However, the cavity function y(r) 
is a much more regular function than the radial distri- 
bution function g{r). In particular, if the interaction 
potential ip{r) diverges in a certain region (as happens 
for HSs with r < a), g(r) vanishes in that region, while 
y{r) keeps being well defined. In that case, y{r) provides 
more physical information than g{r). In this context, it 
is worth recalling the physical interpretation of y{r) as 
being proportional to the probability of accepting the in- 
sertion of an additional test particle at a distance r from 
a given particle with which the test particle is assumed 
not to interact [3ll |. 

Even in the one-dimensional (ID) case, the PS model 
is far from trivial, since interactions are not restricted 
to nearest neighbors and so its exact solution is not 
known [2^ . Useful information, however, can be ob- 
tained in analytical or semi-analytical form in some lim- 
iting cases. Thus, in the combined high-temperature and 
high-density limit the PS model is amenable to an ex- 
act analytical treatment for any number of dimensions 
An interesting consequence of the boundedness of 
the PS potential in the ID case is the plausible existence 
of a fluid-crystal phase transition p6l | , thus providing one 
of the rare examples of phase transitions in ID systems 
[s^ l . In the complementary low-density domain, the cav- 
ity function y{r) has been recently determined to second 
order in density at any temperature [30l |. Comparison 
of the second-order contribution to y{r) with the corre- 
sponding PY and HNC predictions shows that both are 
rather poor inside the core for low and moderate tem- 
peratures (say T* = ksT/e < 1) but the HNC theory 
rapidly improves as the temperature increases. Finally, 
at zero temperature (T* 0), the PS model becomes 
identical with the HS model. Although the correlation 
functions of HSs are not exactly known (except in the ID 
case), the PY theory is analytically solvable for that in- 
teraction model [l| , the corresponding radial distribution 
function exhibitin g a general good agreement with com- 
puter simulations [33|, [SJ]- However, this good agreement 
of the PY g{r) for HSs can be misleading. As mentioned 
before, g(r) vanishes inside the core (r < a) for HSs but 
the cavity function y(r) does not. It turns out that the 
PY y{r) for HSs is much worse for r < a than for r > a 
(even in the ID case [29} ), this deficiency being inherited 
by the PY solution for PSs [30], in which case g(r) ^ 
for r < a. 

The main aim of this work is to construct two sim- 
ple analytical theories for the structural properties of the 
three-dimensional PS fluid. One of them is based on the 
exact solution of the PY equation for HSs, modified in- 
side the core. Since, as noted above, the HS fluid is the 
zero-temperature limit of the PS fluid, this extension of 



the PY solution is expected to be adequate for low and 
moderate temperatures and so we will refer to it as the 
low-temperature (LT) theory. The other approach, here 
referred to as the high-temperature (HT) theory, is a sim- 
ple extension to finite temperatures of the exact asymp- 
totic correlation functions in the limit T* oo. Both 
approximations are compared with our own MC simula- 
tions as well as with numerical solutions of the PY and 
HNC integral equations. It is observed that the LT and 
HT theories complement quite well each other, exhibit- 
ing a good agreement with the simulation data in their 
respective domains of applicability, which are wider than 
what one might have anticipated. It is also found that 
the HNC theory has an excellent behavior, except inside 
the core at low temperatures, while the PY results are 
always very poor inside the core. 

In the next Section some basic properties of the PS 
fluid are presented. Next, the zero-temperature limit is 
considered in Sec. IIIIl where special attention is paid to 
the PY solution for HSs and an improved modified ver- 
sion of it in the overlapping region is proposed. Our LT 
and HT theories are worked out in Sees. lIVl and fVl re- 
spectively. Section IVII shows the comparison with MC 
simulation data and with the PY and HNC results for a 
wide range of representative states. Finally, the paper is 
closed in Sec. IVII] with a summary and discussion of the 
results. 



II. THE PS MODEL. SOME BASIC RESULTS 

Let us consider a three-dimensional fluid made of par- 
ticles interacting via the pairwise PS potential (jl.ip . The 
cavity function y(r) = e^'^^^^^^'^ g{r) is a continuous func- 
tion of r, while the discontinuity of the potential at r = a 
is transferred to the radial distribution function g{r): 

g{r)^y{r)~xy{r)e{l-r), (2.1) 

where 

x = l-e-^/^' (2.2) 

is a temperature-dependent parameter bounded between 
x = (T* oo) and a; = 1 {T* 0), e(z) is the 
Heaviside step function, and henceforth the distance is 
measured in units of cr. The parameter x represents the 
probability of rejecting an overlap of two particles in an 
MC move. Some properties of the Laplace transform G(t) 
of rg(r) are presented in Appendix [Xl 

In the limit of infinitely high temperatures at finite 
densities the height of the potential barrier becomes neg- 
ligible and so the system behaves as an ideal gas. Apart 
from this trivial situation, the exact form of y(r) is not 
known, except (a) to second order in density for any tem- 
perature [30| and (b) in the high-temperature and high- 
density limit [2^ . 
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A. Low-density limit 

To first order in density, the cavity function is exactly 
given by 

y(r) = 1 + ^iix\2 - rf{r + 4)9(2 - r) + 0{r]^), (2.3) 

where 77 = (7r/6)pcr'^, p being the number density. In the 
case of HSs, rj represents the packing fraction and so here 
we will retain that name (although in the PS model there 
is no a priori upper bound for it). From Eq. (12. 3p one 
can easily get 



2/(1) 



l + -77.x2 + 0(r;2), y'(l) 



-rjx 



(2.4) 



y{0) = l + 8r]x'^ + 0{r]^), y' (0) = ~6r]x^ + 0{f]^), (2.5) 

where the prime denotes a derivative with respect to r. 
The second-order contribution in the density expansion 
is also known but it will not be needed here. 



B. High-temperature limit 

In the high-temperature and high-density limit (T* — > 
00, 77 (X), r]x = const), the cavity function has the form 

m 



(2.6) 



y{r) = 1 + xw{r) + Oix'^), 
where w{r) is the inverse Fourier transform of 
967r7ya;(fc cos k — sin fc)^ 



w{k) 



[k^ — 2477a; {k cos k — sin k)] 



(2.7) 



The virial expansion of the function w{r) is worked 
out in Appendix |B1 where it is shown that the series 
converges for 77a; < | . It is also proven in Appendix |B] 
that w(r) has a fourth-order discontinuity at r = 1 and 
a discontinuity of order 2(77, — 1) at r = 77 > 2. 

The denominator in Eq. (|2.7p has a real root if 7;a; > 
770 ~ 1.45. This indicates the existence of a spinodal tran- 
sition (Kirkwood's instability) in the high-temperature 
limit as 77a; — > 770 from below. This instability is pre- 
empted by the freezing transition, which can be esti- 
mated to occur at 7;a; « 0.6 [2^. 



III. ZERO-TEMPERATURE LIMIT 
A. Cavity function inside the core 

In the limit of zero temperature (T* — > 0) the pene- 
trable particles become impenetrable, i.e., they become 
HSs. Therefore, in that limit the radial distribution func- 
tion g(r) vanishes inside the core (r < 1). However, the 
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FIG. 1: (Color online) Plot of the cavity function y{r) inside 
the core (r < 1) for densities pa^ = 0.3, pcr^ — 0.5, and 
pa^ = 0.7. The symbols are MC simulation data for HSs fsij . 
The lines are our MC data for PSs at T* = 0.2 (solid lines), 
T* = 0.13 (dashed-dotted lines), and T* = 0.1 (dashed lines). 



cavity function is still well defined in that region. In con- 
nection with this point, it might be useful to recall two 
zero-separation theorems for HS fluids [1^ : 



In 7/(0) 



1 



knT 



Pcx = Z{p) - 1 + 



" , >Z{p') 
dp - 

P 



1 



-67777(1). 



, (3.1) 



(3.2) 



In Eq. (|3.ip , /icx is the HS excess chemical potential and 
Z = 1 + 4riy{l) is the HS compressibility factor. 

It is interesting to see how the cavity function for PSs 
converges to that of HSs as the temperature decreases. 
Figu re [1] compares our MC simulation data obtained for 
PSs Hi^at T* = 0.2 (1 - a; ~ 7 X lO^^), r* = 0.13 
{1- x ~ 5 X 10-^), and T* = 0.1 (1 - a; ~ 5 x 10"^) 
with the MC results obtained by Labik and Malijevsky 
for HSs |3l| . Our PS data were obtained by measuring 



g{r) directly and then obtaining y{r) inside the core as 
y{r) — g{r)/{l — x). In contrast, the HS data of Ref. [3l| 
were obtained by a particle-insertion method. 

We observe that for a density pa^ — 0.3 a tempera- 
ture T* = 0.2 is low enough to get results practically 
indistinguishable from those of HSs (T* = 0). However, 
as the density increases one needs to keep lowering the 
temperature to get a good collapse. Thus, for pa^ — 0.5 
and pa^ — 0.7 the HS curves are recovered at T* = 0.13 
and T* =0.1, respectively. It is interesting to note that, 
although at T* = 0.2 the probability of accepting an at- 
tempted overlap between two PSs is only 1 — a: ~ 7x10^"^, 
a density p = 0.7 is high enough to make the correspond- 
ing cavity function deviate strongly from that of HSs. 
Figure [T] also shows that at low densities and temper- 
atures (e.g., pa^ = 0.3 and T* = 0.1) the number of 
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overlaps is very small and so the measured y{r) presents 
large statistical fluctuations. 



B. PY solution for HSs 

It is well known that the radial distribution function 
gOr) for HSs is rather well represented by the PY solution 
[III IS [13 • that solution the auxiliary function P{t) 
defined by Eq. (jASp has the explicit rational form 



Pit) = -e- 



l + Lit 



1 + Sit + S2t^ + Sst^ ' 



(3.3) 



where the coefficients can be uni vo cally determined by 
application of Eq. ^M^. This yields [13 [ll] 



Li 



1 + ^?/2 

1 + 27? ' 



(3.4) 



(3.5) 
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FIG. 2; (Color online) Plot of the HS cavity function y{r) 
inside the core (r < 1) for densities pa^ — 0.3 (squares), 
pa^ = 0.5 (triangles), and pa^ = 0.7 (circles). The symbols 
are MC simulation data [3l[|, the solid lines are the PY pre- 
dictions, Eq. (|3.13|l , and the dashed lines are the modified PY 
predictions, Eq. p.l7p . 



^3 



6r]Li - 1 - 27] 



127; 



(3.6) 



(3.7) 



Insertion of Eq. (|3.3p into Eq. (jA5[) gives the Laplace 
transform G{t) for HSs. It can be expressed as 



G(i) = ^F„(t)e 



n=l 



where 

Fn{t) 



t 



1 + Lit 



127? \1 + Sit + S2t^ + S3t\ 
The inverse Laplace transform of Eq. p.Sp yields 



(3.8) 



(3.9) 



^ OO 

air) = - X! ~ ~ n)' 



(3.10) 



n=l 



where /n(r) is the inverse Laplace transform of Fn{t). In 
particular, 

l + 7?/2 , 97?(l+77) 

yPYW = - ypY(l) = - 2^^_^^3 ' (3.11) 



(l-7y)2 



where henceforth we use the label PY in subscripts or 
superscripts to indicate explicit results derived from the 
PY solution for HSs. The contact value ?/py(1) yields the 
virial route to the equation of state and from it one gets 



(3.12) 



Equation p.lOp gives reasonable predictions for the HS 
radial distribution function (33l. [33|. However, the cavity 
function in the region r < 1, which is given by 



ypvir) = 



1 



(1 + 277)2 -677(1 + 77/2)2r 



+ |(1 + 277)V 



(3.13) 



strongly underestimates the correct values, even to sec- 
ond order in density 30]. From Eq. ()3.13p one gets 



ypY(O) 



(1 + 277)2 y^^iO) 



-67/ 



1 + 77/2 



(1 - 77)4 ' ypY(O) V 1 + 277 

(3.14) 

Equations (|3.11[) . (|3.12p . and (|3.14p are not consistent 
with the zero-separation theorems p.ip and (|3.2p . In 
fact, the values of 7/(0) and y'{0) thermodynamically con- 
sistent with 77py(1) are 



ymPY(O) = (1 - 77)2exp 



277(5 - 277) 



2/mPY(0) 



-677 



1 + ^/2 

(1-77)2' 



(3.15) 



(3.16) 



where here the label mPY stands for "modified" PY val- 
ues. One can then modify the PY expression of y{r) for 
< r < 1 as 



2/mPY(r)=7/PY(r)e('-i)'(^°+^^'^), 



(3.17) 



where 77mPY(l) = 2;py(1), 2/mPY(l) = ypY(l)' and the 
coefficients Aq and Ai are determined by imposing con- 
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sistency with Eqs. (|3.15p and (|3.16p . namely 



Ao = In 



ymPY(O) 
2/py(0) ' 



Ai — 2Ao 



ymPY(O) 2/py(0) 



ymPY(O) ypY(O)' 

(3.18) 

Of course, the modified form (|3.17p is not the only 
possible one that can be proposed to satisfy the 
zero-separation theorems with ymPY(l) — ypY(l) and 
y^PY(l) = ypY(l)- The rational behind Eq. (PTZ)) is 
that, since the HS cavity function grows so rapidly inside 
the core (see Fig.[T|), it is convenient to add a polynomial 
correction to Inypy^r) rather than to ypyir). 

Figuredlcompares Eqs. fXT^ and ([XT7)) with MC data 
[3H . It is quite apparent that the PY is very inaccurate in 
the region < r < 1. However, the simple modified form 
(I3.17p . apart from incorporating the values of y{0) and 
y'{0) consistent with the PY contact value t/pY(l) via 
the zero-separation theorems, is quite satisfactory. On 
the other hand, 2/mPY(?') underestimates the true values, 
this effect being especially visible at pa^ — 0.7. This is a 
consequence of the fact that, as is well known, ?/py(1) is 
lower than the true contact value and this is transferred 
to j/mPY(0) through Eq. 

In order to extend Eq. (|3.15p to the PS case in Section 
IIVI it is convenient to express ?/mPY(0) in terms of j/py(1). 
To that end, we first rewrite Eq. p.lSp in the form 



ymPY(O) = yPY(l)e"('')™^-(^\ 



The steps are similar to the ones followed in the ID 
case In a first stage, we keep the relationship (|A5p 
between the Laplace transform G{t) and the auxiliary 
function P{t) and propose an explicit form for the lat- 
ter. In a second stage, we modify the function g{r) in 
the overlapping region in order to satisfy approximate 
extensions of the zero-separation theorems. 



A. First stage 



Our proposal for the function P{t) defined by Eq. (|A5P 



IS 



which reduces to the PY form (|3.3p for HSs if Lo = 1- 
However, the presence of 7^ 1 for T* ^ (i.e., a; 7^ 1) is 
necessary to satisfy the exact physical requirement (jA7p . 
Condition (|A6p allows one to express the coefficients Li, 
Si, 5*2, and 6*3 in terms of Lq: 



1 + 1]/2 



(3.19) 



Si — Li — Lq 



(4.2) 



(4.3) 



where 



1 



l + r,/2 



2(5-277) 



In- 



7] 1 + r]/2 

(3.20) 

It turns out that a{ri) depends very weakly on 77, deviat- 
ing from a(0) = ^ less than 5% for 77 < 0.48. Therefore, 
from a practical point of view we can approximate ()3.15p 
or (jXig)) by 



7/^PY(0)~7;pY(l)e^™--(i). 



(3.21) 



In fact, we have checked that Eqs. (|3.17p and (|3.18p . when 
implemented with Eq. p.2ip instead of Eq. (|3.15p . pro- 
vides curves indistinguishable from the mPY curves plot- 
ted in Fig. [51 except at the density pa^ = 0.7, where Eq. 
p.2ip gives slightly better results. 



IV. THE LOW-TEMPERATURE (LT) 
APPROXIMATION 



S2 — -^Lq — Li, 



(4.4) 



^3 = ^^^^^-l~^^^^o _ (4 5) 

12ri 

Inserting (|4.ip into (|A5p and expanding in powers of 



e % one gets 



G(t) = ^F„(t)e-"*, 



(4.6) 



71=0 



where 



Foit) 



1277 Lo + Sit + S2t^ + Sst^ ' 



(4.7) 



Fn{t) = - 



t {l + Sit + S2t^ + Sst^) {Lo + Lity 



12?/ {Lo + Sit + S2t^ + Sst^r^^ ' 
n>l. (4.8) 



Now we want to generalize to PSs the PY solution for 
HSs described in Sec. lIIIi plus the modification (|3.17p for 
r < 1, so that effects due to finite (but low) temperatures 
are incorporated. The resulting approximation for g{r) 
will be essentially analytical, will reduce to the PY so- 
lution in the limit T* — > 0, and will be better than the 
numerical solution of the genuine PY integral equation 
for PS model, except at high temperatures. 



Therefore, 



00 

9{r) = -Y.fn{r-n)e{r-n), (4.9) 



n=0 



where /n(7') is the inverse Laplace transform of F„ {t). 
Equations (|4.6p and (|4.9p are formally analogous to Eqs. 
and p.lOp . respectively. However, some important 
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differences can be observed: Eqs. (|4.6p and (|4.9p include 
extra contributions associated with n — 0, Eq. (|4.8p dif- 
fers from Eq. ([aj]) . and Eqs. (|4?2|l - (|45| differ from Eqs. 

Obviously, Eqs. g21)-gH) reduce to Eqs. 
([Oll - laJOl) if Lq = 1. Since Lq ^ 1 for PSs, the basic 
difference with respect to the HS case is that now one has 
.9('') = /o('')/^ 7^ for r < 1. More specifically, appli- 
cation of the residue theorem to get the inverse Laplace 
transform of (|4.7p yields 



Ln-1 



t — 1 ^ 



(4.10) 



where Zi {i = 1,2,3) are the three roots of the cubic 
equation Lq + Sit + S2t^ + S-^t^ = 0. In particular. 



Lo-l 



12tj^ ' Si + 2S2Z., + SSsz, 

2—1 ' 



re^% (4.11) 



Lo-l 



12?7^ ' Si + 2S2z^ + 3S3zf 

Z — 1 '' 



(4.12) 



One could also evaluate /o(0) and /g(0) from Eq. (|4.10p . 
but it is easier to do it by a simpler alternative route. 
Taking into account the expansion of Fo{t) for large t, 



Foit) - 
one gets 

Mr) 

Therefore, 

Similarly, 
Flit) = 

so that 

Mr) = 

/i(0) = 



Lo-lt- 



Ss 1277 



S 



2,-1 



+ 0{t-'^), (4.13) 



/o(o) = o, MO) 



Lo-l 
l2vS3 ■ 



(4.14) 



(4.15) 



Li 

^12^ 



S3 12?/ 



+ 0(r'), 

(4.16) 



Lo S2 



S3 



+ 0{r^), (4.17) 



1277^3 



/{(0) = -^^^4^^. (4.18) 



So far, the parameter Lq remains undetermined. We 
now fix it by imposing the continuity of the cavity func- 
tion at r = 1, i.e., j/(l~) = 2/(1^)- From Eq. (|4.9p one 
has 



il-x)-^fo{r), 0<r<l, 



'■^^'^^^ l/o(r)+/i(r-'l). 



1 < r < 2. 



(4.19) 



Therefore, 



l-x 



foil), 2;(l+) = /o(l) + /i(0). (4.20) 



The condition 7/(1 )=7/(l+) yields 

a;/o(l) = (l-a;)/i(0). 



(4.21) 



Insertion of Eqs. and into Eq. ^i^, to- 

gether with Eqs. (|4.2p - (|4.5p . renders a closed transcen- 
dental equation that gives Lq as a function of x and 77. 
This closes the proposal (|4.9p . In the zero-temperature 
limit (x ^ 1), Eq. (I4.2ip becomes /o(l) — 0, whose solu- 
tion is simply Lq — I and then we recover the analytical 
PY solution for HS s. 

While Eq. (|4.2ip satisfies the condition yil^) = y(l^), 
it does not verify y'(l~) = j/'(l+), except at T* = 0. 
In addition, it inherits from the PY solution for HSs a 
poor performance in the overlapping region, even at low 
temperatures. Both deficiencies are remedied by the next 
stage. 



B. Second stage 

In parallel with what we did in Section IIIII in the case 
of the PY solution for HSs [cf. Eq. (|3.17p ]. we propose 
now in the case of PSs an improved version of the ap- 
proximation (|4.9p in which only the radial distribution 
function in the overlapping region is modified, namely 

Q(r)e(l-r) °° 

9ir) = - : ^/„(r-7i)e(r-7i), (4.22) 

n=0 



r 



where Qir) is a third-degree polynomial with the con- 
straint Qil) — to maintain the continuity of ?/(r) at 
r = 1 obtained in the first stage. The degree of Qir) is 
suggested by the structure of the exact cavity function 
to first order in density, Eq. (j2.3p , and by the polynomial 
appearing in the exponential of Eq. p.l7p . Thus, Qir) 
can be written as 

Q{r) = (r - 1) [A + B(r + 2)(r - 1) + Cr{r - 1)] . 

(4.23) 

The parameters A, B, and C will be determined by im- 
posing the continuity condition y'il^) = 7/'(l+), as well 
as approximate expressions for the zero-separation val- 
ues 7/(0) and 7/'(0). According to Eq. (5221), Eq. is 
replaced by 

^" I /o(r) + /i(r-l), l<r<2. ^^-'^^ 
The condition 7/'(l^) = y'(l^) yields 



A = 



/i(0) 



:/o(i) 



(4.25) 



where use has been made of Eq. (|4.2ip and /o(l), /i(0), 
and /{(O) are given by Eqs. (|4.12p and (|4.18p . It remains 
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to determine the coefficients B and C. To tfiat end, note 
that Eq. ((i?^ imphes 



r* r? y(l) Iny(O) Eq. -i/'(0)/y(0) Eq. gST]) 



2/(0) 



l-x 127yS'3 



y'jo) 
yio) 



= A + C-3B 



2^3' 



(4.26) 



(4.27) 



where use has been made of Eq. (|4.14p . Therefore, B and 
C can be expressed in terms of y{0) and y'(0) as foUows: 



(1 - x)Ur]- 



-2/(0) 



C 



y'jo) 

2/(0) 



(4.28) 



-^ + 35+^. (4.29) 



To close the problem we need to propose approximate 
expressions for y{0) and y'(0). Here we take 
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In ?;(0) = — ^ix y{l) + In 
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y{l) + T^{x' ~ x'')y\l) 



y'jo) 

2/(0) 



= —6rix^y{l). 



(4.30) 
(4.31) 



Equations (|4.30p and ()4.3ip are exact to first order in 
density, as can be seen from Eqs. p.4p and (|2.5p . In 
addition, Eqs. ((Oil)) and (|i3T|) reduce to Eqs. ([^^ and 
(|3.2p . respectively, in the zero-temperature limit {x — ^ 1). 
From that point of view, Eqs. (|4.30p and (|4.3ip can be 
viewed as approximate extensions to a; 7^ 1 of the zero- 
separation theorems. A similar extension was proposed 
in the ID case (29| . 

To check the reliability of Eqs. ((13(1)1 and (|OT)) . we 
have tested them against MC simulations for a wide range 
of temperatures and densities. The results are displayed 
in TableHI It can be seen that Eqs. (|T^ and (lOT)) com- 
pare reasonably well with simulations for all the states 
considered. 

This closes the construction of our low-temperature 
(LT) approximation. In summary, for any desired den- 
sity and temperature, the radial distribution function is 
given by Eq. (|4.22p . where fa{r) is provided by Eq. (|4.10p 
and fn{f) with n > 1 are the inverse Laplace transforms 
of the functions (|4.8p . Some needed specific values are 
given by Eqs. (|4.1ip . (|4.12p . and (|4.18p . The parameters 
Li, Si, S2, and S3 are linear functions of Lqi a-s shown by 
Eqs. (|4.2p - (l4.5p . and the latter quantity is the solution of 
Eq. (|4.2ip . Finally, the polynomial Q{r) is given by Eq. 
(|4.23p . where the explicit expressions for the coefficients 
A, B, and C can be found in Eqs. (HSS)), and 

(|4.29p . the two latter supplemented by Eqs. (|4.30p and 
(|4.3ip . where y{l) is given by Eq. (|4.20p . AU the expres- 
sions are analytical, except Eq. ()4.2ip . which must be 
solved numerically to obtain Lq for given 77 and x. Note 



0.1 


0.3 


2.48 


4.85 


5.00 


4.27 


4.46 


0.1 


0.37 


3.16 


7.37 


7.58 


7.30 


7.01 


0.2 


0.3 


2.30 


4.52 


4.57 


5.11 


4.08 


0.3 


0.3 


1.91 


3.53 


3.54 


3.65 


3.20 


0.4 


0.3 


1.67 


2.75 


2.73 


2.94 


2.53 


0.5 


0.3 


1.52 


2.18 


2.14 


2.34 


2.05 


0.5 


0.4 


1.58 


2.85 


2.79 


3.48 


2.84 


0.5 


0.6 


1.59 


3.94 


3.87 


5.66 


4.28 


1.0 


0.4 


1.25 


1.16 


1.10 


1.31 


1.20 


1.0 


0.5 


1.28 


1.38 


1.32 


1.67 


1.53 


1.0 


0.6 


1.29 


1.62 


1.50 


2.31 


1.86 


1.0 


0.8 


1.30 


2.03 


1.84 


2.99 


2.49 


1.7 


0.6 


1.15 


0.72 


0.73 


0.87 


0.82 


1.7 


0.8 


1.16 


0.91 


0.89 


1.20 


1.10 


1.7 


1.2 


1.17 


1.25 


1.16 


1.88 


1.67 


3.0 


0.8 


1.08 


0.37 


0.40 


0.44 


0.42 


3.0 


1.2 


1.09 


0.50 


0.55 


0.73 


0.63 


3.0 


1.8 


1.10 


0.68 


0.73 


0.98 


0.95 


6.0 


3.0 


1.05 


0.31 


0.39 


0.48 


0.45 


20.0 


6.0 


1.01 


0.06 


0.09 


0.07 


0.09 



TABLE I: MC simulation data for y{l), Iny(O), and 
j/'(0)/y(0). The values obtained from the approximations 
(|4.30p and (|4.3H) with the empirical j/(l) are also included. 



that, in order to determine g{r) for r < k — integer, 
where in practice fc = 3 is enough, only the functions 
/„(r) with n < k — 1 are needed. They can be explicitly 
obtained by application of the residue theorem, as illus- 
trated in Eq. ()4.10p for n = 0. Alternatively, one can 
use any of the efficient methods discussed by Abate and 
Whitt [39i| to numerically invert Laplace transforms. 

The LT theory for PSs proposed in this Section is 
within the spirit of the rational-function approximation 
that has been applied in the past to hard-core systems 
[3^1 • In fact, as indicated before, it is straightforward to 
prove that the LT g{r) given by Eq. (I4.22p reduces to the 
FY g{r) in the HS limit {x ^ 1 ^ Lq ^ 1). The HS 
limit of the LT y{r) for r < 1 is presented in Appendix [Cl 
It turns out that the resulting function, Eq. ()C5p . gives 
values that are practically indistinguishable from those 
obtained for HSs by Eqs. ()3.17p and ()3.18p implemented 
with Eq. (|3.2ip . It is also be proven in Appendix [C] that 
our LT approximation ()4.22p for PSs is consistent with 
the exact form ()2.3p of y{r) to first order in density. 



THE HIGH-TEMPERATURE (HT) 
APPROXIMATION 



The exact asymptotic behavior (|2.6p in the limit T* 
00 can be used to propose a simple approximation that 
might be useful for high temperatures, in the same spirit 
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as done in the ID case [29j. The key point of the ap- 
proximation is the assumption that the cavity function 
y(r) depends on the distance r and the state variables x 
and 77 only through the combination xw{r), where w{r) 
is the inverse Fourier transform of Eq. (|2.7p . which de- 
pends on X and rj through the product r]x. This general 
assumption can be formally expressed as 



y{r) = T{xw{r)) , 

where the function J-'{z) remains to be chosen, 
tency with Eq. (|2.6p imposes the constraint 



(5.1) 



Consis- 



(5.2) 



We have found that T{z) = 1 -f ze^ compares with MC 
simulations for high temperatures better than other sim- 
ple possibilities, such as T{z) = 1-1-2;, J-{z) = e^, or 
J-{z) = (1 — z)~^. Therefore, the high-temperature (HT) 
approximation adopted in this paper is 



y(r) = l + 2;u'(r)e^"'(''). 



(5.3) 



We have also checked that Eq. (|5.3p is superior to the 
mean-field approximation c(r) = /(r), where c(r) is the 
direct correlation function and /(r) — —xQ{l — r) is 
the Mayer function. The Ornstein-Zernike relation then 
yields g(r) = 1 + xw{r) + /(r), which implies 



y{r) = 1 + xw{r) 



l-x 



9(1 



(5.4) 



While the mean-field approximation (|5.4p is also consis- 
tent with the exact behavior (|2.6p . it does not belong in 
the class of approximations (|5.ip . In fact, all the approx- 
imations of the form (|5.ip satisfying Eq. (|5.2p are consis- 
tent with the exact low-density behavior (|2.3p , while the 
mean- field approximation (j5.4p is not. 




VI. COMPARISON WITH MC SIMULATIONS 

Here we compare our two heuristic theories with MC 
simulations and with the results obtained by the numer- 
ical solution of two classical integral equation theories: 
the PY and the HNC approximations. 

We first fix the packing fraction at 77 = 0.3 and in- 
crease the temperature from T* —0.1 to T* = 0.5. The 
results for g{r) are shown in Fig. [31 We observe that the 
LT, PY, and HNC theories practically coincide with the 
simulation results outside the core (r > 1), while the HT 
approximation does it for longer distances (r ^2). In any 
case, the interesting region for the PS model is the over- 
lapping one (0 < r < 1) and it is there where the largest 
discrepancies appear. As expected from the behavior of 
y(r) in the zero-t emp erature limit (see Fig. [2]) and also 
for low densities [30|, the PY theory strongly underes- 
timates g{r) in the overlapping region. We will see in 
subsequent figures that this drawback persists for higher 



FIG. 3: (Color online) Plot of the radial distribution function 
at 7? = 0.3 and for temperatures T* = 0.1, T* = 0.2, T* = 0.3, 
T* = 0.4, and T* = 0.5. The circles are MC simulation data, 
while the lines correspond to the LT approximation ( — ), the 

HT approximation ( ), the PY theory ( ), and the HNC 

theory (- • -). 



densities and/or temperatures. As for the HNC theory, 
it tends to overestimate g{r) for r < 1, this effect being 
dramatic at very low temperatures (say T* = 0.1 and 
T* — 0.2) but rapidly diminishing its importance as the 
temperature increases. In fact, the HNC theory performs 
a remarkable good job at T* = 0.4 and T* = 0.5, except 
near the origin. Concerning the HT approximation, it is 
not surprising that it is so poor, given the range of tem- 
peratures considered in Fig. [3l The most noticeable fea- 
ture is that our LT approximation, which only requires 
to solve numerically Eq. (I4.2ip . exhibits a surprisingly 
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MC 
-LT 
-HT 
■PY 
•HNC 



1.5 



0.1 



FIG. 4: (Color online) Plot of the radial distribution function 
at T* = 0.5 and for densities ri = 0.3, 77 = 0.4, and r] = 0.6. 
The circles are MC simulation data, while the lines correspond 

to the LT approximation ( — ), the HT approximation ( ), 

the PY theory ( ), and the HNC theory (- • -). 
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good agreement with simulation for the states consid- 
ered. Therefore, it achieves its goal of extending to low 
temperatures the good role played in the case of HSs by 
the (modified) PY theory described in Section HTll While 
the temperature T* = 0.1 (1 — x ~ 5 x 10~^) is so low 
that the PS fluid with a packing fraction 77 — 0.3 behaves 
practically as the HS system, and so ^(r) w for r < 1, 
for temperatures larger than T* = 0.2 (1 — a; « 7 x 10^^) 
the number of overlapped pairs becomes observable. The 
state 2 = 0.3 and T* = 0.2 was already considered in 
Refs. [1^ and [2l|, where it was also found that the PY 
and HNC theories were highly insufficient in that case. 

In Figs, in m m and [7| we consider the temperatures 
T* = 0.5, T* ^ l,T* = 1.7, and T* = 3, respectively, 
and several increasing densities in each case. The val- 
ues of the densities are chosen as to make the product rjx 
(which is the relevant parameter in the high-temperature 
limit ^26]) ranging from -qx ~ 0.25 to rjx ~ 0.5. The 17 in- 
dependent states considered in Figs. [SHZ] are represented 
in Fig. [H 

It is observed in Figs. SHT] that the HNC theory pro- 
vides excellent predictions for T* > 0.5 at any density. 
In contrast, as anticipated before, the PY theory keeps 
being very limited at short distances, although it tends 
to improve as the temperature increases. With respect 
to the LT approach, Figs. [3H7] confirm that it starts fail- 
ing as the temperature grows. The density also plays an 



FIG. 5: (Color online) Plot of the radial distribution function 
at r* = 1 and for densities rj = 0.4, 77 = 0.5, 77 = 0.6, and 
r] = 0.8. The circles are MC simulation data, while the lines 
correspond to the LT approximation ( — ), the HT approxi- 
mation ( ), the PY theory ( ), and the HNC theory 

(- ■ -)• 



important role since, at a given temperature, the quality 
of the LT predictions worsens with increasing density. In 
any case, it is remarkable that the LT theory is rather 
good for temperatures as high as T* — 1.7, provided 
that the density is not large (say 77 ^ 0.6). However, the 
LT theory visibly deteriorates for higher temperatures 
and densities, developing an artificial local minimum at a 
distance r < 1. In contrast, the HT theory, as expected, 
becomes more accurate when the temperature and the 
density increase. In fact at T* = 3 the HT curves prac- 
tically coincide with the HNC ones. 

As shown in Figs. [3HZl the LT and HT theories com- 
plement each other: the former becomes more accurate 
as the temperature and/or the density decrease, while 
the opposite happens for the HT theory. In fact, both 
heuristic theories meet and become practically equiva- 
lent at certain intermediate temperatures and densities, 
as illustrated by the top panel of Fig.[Sl In order to char- 
acterize the "basins of attraction" of each theory in the 
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1.25 




FIG. 6: (Color online) Plot of the radial distribution function 
at T* = 1.7 and for densities 77 = 0.6, r) — 0.8, and rj = 1.2. 
The circles are MC simulation data, while the lines correspond 

to the LT approximation ( — ), the HT approximation ( ), 

the PY theory ( ), and the HNC theory (- • -). 



density-temperature plane, let us define the locus of state 
points where the contact quantity y(l) = g{l~^) takes the 
same value in both approximations. This locus, which is 
plotted in Fig. [51 separates the basins of each approxima- 
tion. As one departs from this curve, the quality of the 
corresponding approximation (either LT or HT) signifi- 
cantly improves, as illustrated by Fig. [3] in the LT case 
and by Fig. [7] in the HT case. Figure [5] also includes the 
Kirkwood's instability line rjx = 1.45, above which the 
function w{r) is not defined. In the high-temperature 
limit the freezing transition has been estimated to occur 
at r/x ~ 0.6 [2^. The extrapolation of this value to any 
temperature provides an (admittedly rough) estimate of 
the freezing curve (dashed-dotted line in Fig. [5]) . There- 
fore the fluid is stable approximately below that curve 
and it is in that region where the domains of the LT and 
HT approximations are meaningful. 



VII. SUMMARY AND CONCLUDING 
REMARKS 

In this paper we have investigated the structural prop- 
erties of the PS model by computational, numerical, and 
analytical tools. In the limit of asymptotically high tem- 
peratures the linear chain diagrams contributing to the 
virial expansion of the cavity function dominate, their 



1.25 




FIG. 7: (Color online) Plot of the radial distribution function 
at T* = 3 and for densities 77 = 0.8, rj = 1.2, and rj = 1.8. The 
circles are MC simulation data, while the lines correspond to 

the LT approximation ( — ), the HT approximation ( ), 

the PY theory ( ), and the HNC theory (- ■ -). 




FIG. 8: (Color online) The circles represent the states con- 
sidered in Figs. [SHZl The thick solid line represents the locus 
of points where the LT and HT approximations predict the 
same value of the contact quantity y{l) = g{l^)- The LT ap- 
proximation is more accurate in the shaded region to the left 
of the curve, while the HT approximation is more accurate 
to the right of the curve. The dashed-dotted line is a rough 
estimate of the freezing curve, while the dashed line repre- 
sents the threshold curve above which the HT theory ceases 
to exist. 



11 



sum giving rise to Eqs. (|2.6p and (|2.7p . The auxiliary 
function w{r) depends on density and temperature only 
through the scaled quantity rjx, where x is defined by 
Eq. (|2.2p . and its series expansion converges uniformly 
for rjx < ^. 

In the opposite zero-temperature limit, the system be- 
comes an HS fluid, whose expression for g(r) in the PY 
approximation is known exactly in Laplace space. Since 
the particles are impenetrable in this limit, g(r) vanishes 
for r < 1. However, the cavity function ?;(7') remains 
different from zero, its determination being important to 
understand the behavior of g{r) for low (but non-zero) 
temperatures, as shown in Fig. [T] We have seen that the 
cavity function for r < 1 predicted by the PY theory 
for HSs is extremely poor (cf. Fig. [2]), this effect being a 
precursor of the tendency of the PY theory to dramati- 
cally underestimate the penetrability phenomenon at fi- 
nite temperatures. However, the modified version (j3.17p . 
implemented with the consistent zero-separation values 
p.lSp and (j3.16p . represents a significant improvement. 
Slightly better results are obtained when Eq. p.lSp is 
replaced by Eq. ([X^ . 

Based on the explicit expression of g{r) obtained from 
the PY approximation in the HS limit, complemented 
by Eq. (j3.17p . we have proposed an approximate theory 
which, by construction, should be accurate for very low 
temperatures (LT). This LT approximation is given by 
Eq. (|4.22p . where the Laplace transforms of the func- 
tions fn{r) can be found in Eqs. (|4.7p and (|4.8p . and 
the polynomial Q{r) is given by Eq. (|4.23p . The coef- 
ficients appearing in those functions are determined by 
imposing consistency conditions, namely the physical re- 
quirement (|A3p . the continuity ofy{r) and y'{r) at r = 1, 
and prescribed expressions for the zero-separation values 
2/(0) and y'{Q)- For the latter two quantities we have 
proposed approximate extensions to finite temperatures 
of the zero-separation theorems holding in the case of 
HSs. Those extensions, which are given by Eqs. (|4.30p 
and (|4.3ip . are reasonably accurate for wide ranges of 
temperature and density, as shown in Table HI All the co- 
efficients are explicitly expressed in terms of one of them 
(chosen to be Lq), which is obtained as the solution of 
a closed transcendental equation, Eq. (|4.2ip . stemming 
from the continuity condition y{l^) — y{l^). 

The LT theory has been complemented by a sim- 
pler approximation constructed by exploiting the ex- 
act asymptotic behavior of the correlation functions for 
high temperatures. Among several possibilities, we have 
taken Eq. (|5.3p as our high-temperature (HT) approxi- 
mation, which agrees with simulation results better than 
the mean-field approximation (|5.4p . 

In order to test the LT and HT theories from moder- 
ately low to moderately high temperatures, we have car- 
ried out MC simulations for the 17 states represented in 
Fig. [H We have also numerically solved the PY and HNC 
integral equations for each one of those states. The re- 
sults show that the LT theory is the best one for T* < 0.4 
and does a very good job for higher temperatures, pro- 



vided the density is not too large. It is a reasonable 
approximation even at a temperature T* = 1.7 and a 
density ?/ = 0.6. However, it rapidly deteriorates for 
higher temperatures (for instance, at T* = 3). As the 
LT approximation worsens, the HT theory improves and 
becomes very accurate for T* > 1.7. Figure [5] shows the 
regions in the parameter space rj-T* where each theory 
is expected to be more reliable 

Regarding the classical integral equation theories, we 
have found that the HNC approximation is excellent, ex- 
cept for low temperatures (T* < 0.4), its performance 
increasing with temperature. This is similar to the situ- 
ation in ID [2^ and agrees with recent results obtained 
for the generalized exponential model with n = 4 . On 
the other hand, the PY theory is generally very poor near 
r = 0, especially for low temperatures. 

The importance of having analytical treatments for the 
structural properties of fluids at one's disposal cannot be 
overemphasized. From that point of view, we believe 
that the LT and HT theories constructed here can be 
useful to describe the spatial correlation functions of the 
PS model without the burden of numerically solving in- 
tegral equations. Both theories complement each other 
and thus their combined use covers satisfactorily well the 
whole range of densities and temperatures. Of course, it 
would be desirable to have a unique analytical approxi- 
mation equally accurate for both low and high tempera- 
tures. However, this does not seem to be an easy task at 
all. As a matter of fact, the PS model represents an inter- 
esting benchmark for liquid theory since it encompasses 
different regimes: the HS system at very low tempera- 
tures and finite densities, the ideal gas at very high tem- 
peratures and finite densities, and the mean-field system 
at very high temperatures and densities. 
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APPENDIX A: LAPLACE TRANSFORM 

The Laplace transform of rg(r) is 

/■OO 

G{t) = / dre^'^Vg(r). (Al) 
Jq 
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The relationship between G{t) and the Laplace trans- 
form H{t) of rh{r), where h(r) = g{r) — 1 is the total 
correlation function, is 



Git) = ^ + H{t). 



(A2) 



The isothermal compressibility is directly related to 
dH {t) / dt\^^Q and is finite. Therefore, H{t) must be finite 
at t = 0, so that the small-t behavior of G{t) must be 



G(i)-^[l + 0(t') 



(A3) 



The large-t behavior of G{t) is also easy to find. 
Near the origin, Eq. (|2.ip implies that g{r) = (1 — 
x) [y(0) + y'(0)r + 0(r2)] , so that 



1 -X 



= ^ [y(0) + 22/'(0)t-i + Oit-^)] . (A4) 

Now we define an auxiliary function P{t) by 
t Pit) 



Thus, Eqs. (Mil and jXI)) are equivalent to 
P(<) = -l-^fVO(t^), 

1 - X 



(A5) 



(A6) 



P{t) = 127?^ [2/(0) + 2j/'(0)t-i + 0(t-^)] , (A7) 



respectively. 



with 



Wn{r) 



24" 
127rr 



dk k sin(A;r) 



k cos k — sin fc \ " 
fc3 



(B5) 

Application of the residue theorem gives, after some stan- 
dard but lengthy manipulations, the result 

, [(»-l)/2] 

Wnir)^- J2 (^-2™-r-)2(""^¥,l™Hr)e(n-2m-r). 



Tn=0 



(B6) 

Here, [k] denotes the integer part of k and V^'' (r) is the 
polynomial of degree n given by 



p=0 



(B7) 



with the coefficients 



12"-in! /3n - 2\ 



(3n - 2)! V P / ^ \q. 



(m — p + g) ! (ri — TO — g) ! 
From Eq. (|B6p it is possible to check the relationship 
w„(l) w„+i(0) 



(B8) 



(n+1) 



= 0, 71 > 2, 



(B9) 



which allows one to prove the thermodynamic consis- 
tency between the virial and energy routes to the equa- 
tion of state [1^ . The first few functions Wn (r) are 

u)2(r) = i(r-2)2(r + 4)9(2 -r), (BIO) 



APPENDIX B: VIRIAL EXPANSION OF THE 
FUNCTION w{r) 

The inverse Fourier transform of w{k) is 



w{r) 



277^7 



dkksm{kr)w{k). (Bl) 



Note now that Eq. (12.71) can be rewritten as 



6rix 1 _ F{k) 6??a; 
where we have called 

k cos k — sin k 



= ^^E^"(^)' (B2) 



F{k) = 24r/a;- 



/c3 



Therefore, Eq. (jBlj) becomes 



(B3) 



(B4) 



3(r-l)\„3 



35r 
35r 



(r^ + 4r^ - 53r - 162)6(1 - r) 
(r^ + 12r2 + 27r- 6)9(3 -r). 



(Bll) 



W4{r) 



{r^2f 
525r 



(r* -f 12r^ - 96r2 - 1232r - 2304) 



x9(2-r) 



(r - 4)^ 



2100r 

-H224r- 144)6(4 -r) 



{r^ + 24r^ + 156r^ 



(B12) 



The series (jB4[) converges for r]x < ^ because w{k) 
is singular at fc = when -qx takes the negative value 
rjx — — i. By an adequate reordering of terms, Eq. (|B4p 
can be rewritten as 

oo 

w{r) = -6(1 - r) + r"^ E - 

n=l 

oo 

xEMn-r)^]"^'"-'pi:LM(B13) 



n=2 
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Therefore, w(r) has a fourth-order discontinuity at r = 1 
and a discontinuity of order 2{n — 1) at r — n > 2. For 
instance, 

w""il+) - w""il-) = 432(^7a;)^ (B14) 



where the coefficients A, B, and C of the polynomial 
Q{r) reduce to [cf. Eqs. Km . gSS]), (l4:28ll - (l43T|) ] 



2/py(1) /o(1)' 



(C6) 



w"(2+) -w"{2-) = -677a;, 



w""{3+)-w""{3-)=A8{r^x)' 



(B15) 



(B16) 



1 , 1 



■In 



ymPY(0)(l-r?)^ 
(l + 2?7)is: 



(C7) 



APPENDIX C: SOME CONSEQUENCES OF THE 
LT THEORY 



1. Zero-temperature limit of the cavity function 
inside the core 



According to Eq. (|4.24p . the cavity function inside the 
core is 



y(r) = (l-x)-i^e«M. 



(CI) 



ymPY(O) " ■ ' 1 - ?/ 

2. Low-density limit 



(C8) 



To first order in density, the parameter Lq has the form 
Lq = Loo+LoiV+^iv^)- The coefficients Lqo and Lqi are 
determined by expanding Fo(i) and Fi{t) to first order 
in 77 and imposing the condition y{l~) ~ ?/(l^) to that 
order. This yields, after some algebra. 



Now we will take the limit e = 1 — a; = e 0. Let 

us first write Lq = Ke + 0{e^). Thus, Eq. (|iTUl) 
becomes fo{r) Kefo{r), where 



127/ ^ 5i + 252 + iS^zj ' 



(C2) 



the parameters Li, ^i, 5*2, and being now given by 
Eqs. (I3.4|) - (I3.7|) . Near the origin, following steps similar 
to those of Eq. (|4.14p , one obtains 



Ur) 



1 + 277 

(1-77)2 



3 7/ 

1-77 (1 - 77) 



3r/2 2 
r 



(C3) 

According to Eq. (|4.2ip . the coefficient K is explicitly 
given by 



ypY(l) 

Therefore, Eq. (jCip becomes 



(C4) 



(C5) 



Lo = 2;+i77a;(l-a;)(3-5a;) + C'(r/2). (C9) 



The functions /o(7') and /i(r) are given by 



/o(f) = (1 - x)r 



1 + -77a;(5 + 5a; - 6r + r^) + 0(77^) 



(CIO) 



/o(r-) + /i(7^-l) = 7^ 



l + -7/a;2(r-2)2(r + 4) + C>(772) 



Finally, use of Eqs. (g^Sl), g^S]), and ([0^1 gives 



(Cll) 



A=-77a;(l-3a:) + 0(772), 



(C12) 



B^--r^x{\~x) + 0{if), C = 0(772) (C13) 



It is then straightforward to check that Eq. (|4.22|) reduces 
to Eq. 
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